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We study the dynamics of rigid spheres embedded in viscoelastic media and address two questions 
of importance to microrheology. First we calculate the complete response to an external force of a 
single bead in a homogeneous elastic network viscously coupled to an incompressible fluid. From 
this response function we find the frequency range where the standard assumptions of microrheology 
are valid. Second we study fiuctuations when embedded spheres perturb the media around them 
and show that mutual fiuctuations of two separated spheres provide a more accurate determination 
of the complex shear modulus than do the fluctuations of a single sphere. 
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Microrheology is an important experimental probe of 
the viscoelastic properties of soft materials Unlike 
more traditional macrorheology, in which a sample is 
subjected to an externally imposed uniform shear strain, 
microrheology relies on the Brownian fluctuations of the 
micron-sized beads dispersed in the sample to assess the 
viscoelastic response function (complex shear modulus), 
G{u>). The principal advantages of this technique are that 
it can be used for the detailed study of materials that 
cannot be produced in bulk quantities and that it can 
be used to probe the local properties of rheologically in- 
homogeneous materials. Because of these two strengths, 
microrheology promises to open a new window on cel- 
lular biology by facilitating the study of the rheological 
properties of intra-cellular structures in living cells. In 
addition, this technique is currently being used to study 
various soft biomaterials 

In a typical microrheology experiment, the time- 
dependent position correlation function of individual 
probe particles is measured either by light scattering 
1^ or by direct real-space imaging This cor- 

relation function provides a complete description (via 
the Fluctuation-Dissipation Theorem) of the frequency- 
dependent response of the probe particles to an external 
force. If inertial effects are ignored, the rheological prop- 
erties of a Newtonian fluid are completely determined by 
a single quantity, its viscosity rj. The displacement of 
a spherical particle, u(cj), of radius a in response to a 
force f(w) at frequency uj in such a fluid is given by the 
standard Stokes-Einstein relation: 
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where G{uj) = —iujrj is the complex shear modulus. A 
natural hypothesis is that this relation can be gener- 
alized to rheologically complex materials in which G{lu), 
the complex shear modulus has both storage (real) and 
loss (imaginary) components. We will refer to this ex- 
tension of Eq. (|l|) as the Generalized Stokes-Einstein re- 
lation (GSER). 

In this Letter we address two basic questions regarding 
the interpretation of microrheological data: (1) In a ho- 



mogeneous viscoelastic medium does the GSER provide 
the correct response function to an applied force? and, 
(2) if the introduction of the probe particles perturbs the 
local rheological properties of the medium, how does one 
extract the unperturbed, bulk rheological properties from 
the data? Recent experiments suggest that, at least in 
certain systems, a discrepancy exits between the macro- 
and microrheological measurement of the shear modu- 
lus 1^ making this question one of current interest. 

To address the first of these questions, we use a model 
viscoelastic medium consisting of an elastic network that 
is viscously coupled to a fluid in which the network is 
embedded In the second half of this letter, we ap- 

proach the problem of rheological inhomogeneities and 
explicitly show that inter-particle position correlations 
are insensitive to the local particle environment, and 
therefore, provide a more reliable probe of the properties 
of the bulk material than do single particle fluctuations 
as proposed in reference 

Our model viscoelastic medium consists of an elastic 
network, characterized by a displacement variable u, that 
is viscously coupled via a friction coefficient T to an in- 
compressible Newtonian fiuid characterized by a velocity 
field V. In the absence of viscous coupling, u obeys the 
standard equation for an isotropic, elastic, compressible 
medium with Lame coefficients A, /i, and v obeys the 
incompressible Navier-Stokes equation with viscosity rj. 
With friction included, the equations for u and v [Q are 

pii - /LtV^u - (A + /i) V (V • u) = -r (u - v) + (2) 
PFV - TyV^v VP = r (u - v) fv (3) 
V • V = 0, (4) 

where P is the pressure and fu and fv are, respectively, 
the force densities exerted on u and v by the embedded 
beads. The friction coefficient F is estimated by consider- 
ing a uniform displacement of the network relative to the 
fluid at constant relative velocity v. The friction force 
per unit volume, Fv, is equal to the friction force rj^v 
on a strand of the network of length equal to the mesh 
size ^ divided by the volume per strand. The result 
is F ~ 7]/^^. 
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Our goal is to calculate the frequency-dependent dis- 
placement compliance a{Lu) relating bead displacement 
r(w) to the external force F{ui) imposed on it: 



r{uj) = a{Lu)F{LLi), 



(5) 



and to determine under which conditions, if any, the 
GSER, a{Lu) = 1/ (diTaG^io)) applies, i.e. under what 
conditions measurements of the displacement of an indi- 
vidual particle provide a direct measure of the complex 
shear modulus of the two-fluid medium. The complete so- 
lution to this problem requires solving Eqs. with time 
derivatives replaced by —ito, fu and fv equal to zero, and 
with boundary conditions that u{uj) = v/{—iuj) = r((jj) 
at the surface of the sphere. The resulting functions 
u(x, w) and v(x, w) can then be used to calculate the 
stress at the surface of the bead and by integration the to- 
tal force Fb(w) exerted on the medium by the bead. New- 
ton's equation for a bead of mass M, —uj'^Mu + F\^{uj) = 
F{uj) then determines a{uj). This procedure is laborious 
at best, and we will apply a slightly less rigorous one. 
We localize the bead-medium forces f^, {a — u, v) on 
the bead by setting w) = F„(a;)6(|k| — fcmax) where 
fa(k, ui) is the Fourier transform of fQ(x, t), Fa{uj) is the 
integrated force exerted by the bead, fcmax — 7r/2a, and 
Q{x) is the unit step function. The total force exerted 
on the bead by the medium is Fb(ijj) = Fu(w) -|- Fv(w) 
and Newton's equation for a bead is the same as above. 

Our procedure is to use Eqs. Ij-^ to calculate u(k, lu) 
and v(k, w) in terms of fu(k, w) and fv(k, w), and then 
to calculate by integration over k, the network displace- 
ment r(ijj) and fluid velocity w(w) at the bead in terms of 
Fa{io). We then require that the bead, the network, and 
the fluid all move together at the bead, i.e. that r{uj) be 
the bead displacement and w(w) = —iojr{u}) its velocity. 
This constraint on r(a;),w(a;) imposes a particular ratio 
between Fu(a;) and Fv(w) that allows us to obtain a lin- 
ear relation between r{uj) and Fb(w). When applied to 
a sphere in a Newtonian fluid and expanded in powers 
of —iuj, this procedure reproduces correctly the constant 
and \/—iu! contribution to a~^{uj) / {—iuj) and the —itu in- 
ertial contribution with a slightly different prefactor. We 
expect similar accuracy for the current problem. Our 
result for a{uj) can be expressed as 
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where we have introduced the complex shear modulus of 
the material: G{uj) — ^—iujr] and the cross-over function, 
H defined by 



H{x) = 1 - / dz [1 



(7) 



as well as the frequency scale: lob = (2/i + A) /(a^F). In 
this result we assume that the mass density of the elas- 
tic network is significantly lower than that of the fluid. 



p/pf <^ 1, owing to the open structure of the network. 
Consequently, in Eq. we have set p = 0. We have 
also introduced the functions J and X which we discuss 
briefly below. A more complete analysis of this result will 
be published elsewhere |6|. In order for the result given 
by Eq. to reduce to the GSER we must find that (at 
least for some frequency range) 7? « 0, J « -^X, and 
/3b(w) = (2pba^w^) /9G{uj) <^ 1 where pb is the mass 
density of the bead. 

First we consider H. From Eq. we note that 
H{x) goes to zero as 1/x for x ^ 1, so for frequencies 
large compared to wb, we can neglect this term. From 
an examination of the hydrodynamic modes of the sys- 
tem, the physical interpretation of this result is clear. 
The frequency scale ujb is the decay time for the over- 
damped longitudinal compression mode of the system at 
the length scale of the bead. In this mode the network 
undergoes a compressional wave while the fluid drains 
from the denser parts of the network to the more rarefied 
parts. The H function, therefore, represents a correction 
to the microrheological measurements due to the exci- 
tation of longitudinal degrees of freedom in the system. 
Whereas in the macrorheological experiment the applied 
strain is pure shear, in the microrheological experiment 
the probe particle responds to all the thermally excited 
modes of the system including the longitudinal compres- 
sion modes of the elastic network. At frequencies higher 
than lub, however, the network "locks in" with incom- 
pressible fluid thereby eliminating the former's longitu- 
dinal modes and bringing the microrheological measure- 
ment into closer correspondence with standard rheology. 
The elimination of the so-called free-draining (longitudi- 
nal) mode at large uj has been discussed previously Q . 

We now consider the function J{uj). Its form is 
controlled by two dimensionless parameters: Pf{^) — 
iuj'^PFa'^ / [G{lj)'k^] and S = (^/a) . The parameter /3f is 
formed by the square of the ratio of the sphere's radius to 
the inertial decay length in the medium and measures 
the importance of fluid inertial effects in the compliance. 
The second parameter, S, simply measures the ratio of 
the network mesh size to the sphere radius. In the limit 
that both /3f, <5 <C 1 the function J{uj) reduces to —X{uj). 
Since pf ~ Pb, Ph and /3f are of the same order and both 
wiU be small for u < lo* with /3b (u;*) ~ /3f(w*) ~ 1. Our 
approximate calculation is expected to reproduce the ex- 
act result for w < w*, so our estimate of the region of 
validity of the GSER should be correct. 

In typical experiments |^ , the probe sphere is taken to 
be orders of magnitude larger than the mesh size so we 
may safely assume that (5^1. For experiments on actin 
Q with a sphere size of 1 pm, (3f remains small up to fre- 
quencies on the order of SOkHz. A similar estimation of 
the lock-in frequency yields wb ^ lOHz. Thus in typical 
experiments there remains a significant frequency win- 
dow, ljb < lu < lu* , where the response function of the 
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probe particle to an applied force is well approximated 
by the GSER. This model calculation reveals the range 
of validity of the GSER for typical experiments on soft 
materials; furthermore it presents a quantitive prediction 
of the form of the compliance in frequency regimes where 
the GSER does not hold §. 

We now turn to the issue of rheological heterogeni- 
ties introduced by the beads themselves. We imagine 
a medium characterized by a homogeneous frequency- 
dependent elastic-constant tensor. The introduction of 
spherical probe particles perturbs the medium in the 
vicinity of these particles and leads to a spatially inho- 
mogeneous elastic constant tensor Kijki{:>c,Lu). Assum- 
ing that the stress-strain relation remains local, that the 
frequency regime {lob < lu < lu* for our coupled net- 
work) is such that the medium can be characterized by a 
single, frequency-dependent elastic constant tensor, and 
that inertial terms can be neglected, the equation for the 
displacement variables is 

- dj {Kijki{:x.,uj)dkUi) = /i(x,a;), (8) 

where fi (x, uj) is the force density that acts on the sur- 
face of the particles. The displacement responses of the 
collection of particles to forces upon them can described 
by a compliance tensor 

i?r(c.') = a[;'")(^)^;™(^), (9) 

where Rf is the displacement vector of the n^^ particle 
and F™ the force on the m}^ particle. We ask which com- 
ponents of the compliance tensor depend on the bead- 
imposed inhomogeneities of Kijki (x, to) and which, if any, 
depend only on the bulk homogeneous part? 

To answer this question, it is useful to consider first 
the simpler but related problem of determining the bulk 
dielectric constant of a medium by measuring the self 
and mutual capacitances of metal spheres whose presence 
perturbs the dielectric constant in their vicinity. If the 
dielectric constant e(x, u) remains local and frequencies 
are such that transverse electric fields can be ignored, 
then the potential 0(x, lo) satisfies 

- V • (e(x, w)V0(x, Lo)) = 47rp(x), (10) 

where /o(x) is the charge density at x. It is clear from 
Eqs. (||) and (|l^) that there is an analogy between the 
electrical and rheological problems with the identifica- 
tion: (f) < — > u, e < — > Kijki, and p < — > f. The total 
charge Q on a metal sphere is the analog of the total 
force F on a bead in the viscoelastic medium. The in- 
verse capacitance tensor defined by 

4'n — C'nmQm, (11) 

where (/>„ is the potential on bead n and Qm is the total 
charge on bead m, is the analog of the compliance tensor. 



To keep our calculation simple, we consider two con- 
ducting spheres of radius a separated by a distance r in 
a medium of dielectric constant e. Each sphere perturbs 
the medium locally, producing a spherical region of ra- 
dius a' with a dielectric constant e as shown in Fig. (|^). 
To leading order in a/r, the inverse self-capacitance is 

This result shows that fluctuations of a single bead are 
sensitive to the local environment around the bead and 
that, therefore, they do not measure directly the bulk 
dielectric constant, e. The inverse mutual capacitance, 

C.V^J^(l + ^'(^)), (13) 

depends, however, only on the bulk dielectric constant 
to leading order in a/r. Thus correlated voltage fluctua- 
tions (w)(/)2(— ti))) — 2 {T/lo) ImC]^^(w) yield a direct 
measurement of e(aj) provided the beads are far enough 
apart that proportional to 1/r. 

Given the formal analogy between the electric and me- 
chanical problems, it is reasonable to assume that dis- 
placement fluctuations of a single bead do not provide a 
direct measure of the bulk rheological properties whereas 
correlated fluctuations of two beads do, provided the 
beads are far enough apart. Indeed, we will find this to 
be the case, however, the vectorial nature of the elastic 
problem leads to some complications. 

We begin by considering a single sphere of radius a that 
perturbs the elastic medium in which it is embedded out 
to a radius a' . For r > a' , the medium is characterized 
by bulk Lame coefficients /i(a;) and X{uj). For r < a' the 
Lame coefficients are /2(w) and A(w). Elastic displace- 
ments Uinncr in the inner (r < a') and Uoutcr h^ the outer 
(r > a') regions satisfy the equation 

^V2u + AV(V-u) = 0, (14) 

where the Lame constants take their appropriate values 
in each region. Putting the force applied to the sphere 
in the z direction, the most general solution for u is 

a A a^B 
u ~ — {vf cos 9 + z) -\ — (3r cos 9 — z) + (15) 

Dr^ 

+Cz -\ — (err cos 9 — z) 

where A, B, C, and D are constants. The constants 
and a depend only on the local Lame constants. The so- 
lution for the strain field as written in Eq. ( p^ includes a 
superposition of a part that decays with distance as 1/r 
and a dipolar term. These first two terms are accompa- 
nied by two other solutions: a constant shift and a term 
growing with distance from the sphere. The later two 
terms cannot occur in Uoutor as that field must go to zero 



3 



at large distances from the sphere. Thus in the bulk so- 
lution for the strain we have two undetermined constants 
{A and B) while the inner solution (uinner) has four un- 
determined constants. The application of the boundary 
conditions at infinity has reduced the problem to find- 
ing six constants. Due to the rigidity of the sphere, the 
displacement is fixed at its surface. This boundary con- 
dition contributes two more constraints. The remaining 
four conditions come from strain field continuity at the 
interface of the two elastic media (r = a') and the conti- 
nuity of the two components of the stress tensor at that 
interface, arr and are- The problem is now completely 
determined. 
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FIG. 1. I. Schematic of a system in which conducting 
spheres of radius o embedded in a medium of dielectric con- 
stant e are surrounded by concentric spheres of radius a' with 
dielectric constant e. II. Similar schematic of rigid spheres 
embedded in an elastic medium with Lame coefficients /i, A 
surrounded by concentric spheres with Lame coefficients p., X. 

The complete solution shows that the self-component, 
a}j , of the compliance tensor depends in a complex way 
on the local Lame coefficients fl and A. Thus as in the 
electrical case, fluctuations of a single bead will not yield 
reliable measurements of bulk rheology unless a' /a — 1 <C 
1 or /i, A do not differ significantly from /i, A. 

To compute the cross component, afj{uj), of the com- 
pliance tensor relating displacements of bead 2 to forces 
on bead 1, we observe that bead 2 will follow the dis- 
placement field produced by bead 1 at separations r large 
compared to a. Thus afj is simply the coefficient of 
in the displacement field of bead 1. At large r, only the 
first term in Eq. (|l5|) survives. The coefficient A in this 
term is determined by a global property of the stress field 
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where the integral is over any closed surface surrounding 
the sphere. Only the l/r part of the displacement field 
contributes to this integral. From this constraint we can 
calculate Aoutcr, the coefficient of the the first term in 
Eq. ( p^ ) in the outer region (r > a'). This coefficient 



is linear in F^. From this, we find that the compliance 
tensor can be decomposed into parts, a\\, parallel to the 
vector r separating the two beads and, q;j_, perpendicular 

to r: afj(w) = a\\fifj -t- a±_ {Sij — fifj) with 
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a± 



47rr/i(a;) 

1 X{uj) 



3/i(w) 



(17) 



(18) 



87rr/i(ci;) A(a;) + 2/i(a;) 

Thus fluctuations parallel to the separation vector de- 
pend only on the shear modulus, /J.(w) = G(cij), whereas 
those perpendicular to the line of centers depend on both 
A and /i. In the incompressible limit, a±/a\^ = 1/2, 
which is identical to the ratio of the parallel and per- 
pendicular diffusivities of two spheres with (incompress- 
ible) hydrodynamic interactions |8|, in agreement with 
recent experimental results on two-point microrheology 
in a viscous liquid The experimental determination 
of this ratio in viscoelastic materials can be used to test 
for compressibility effects at the frequencies relevant to 
the experiment. 

The combination of single-particle and two-particle po- 
sition correlations provide data about both the local en- 
vironment of the probe particle and the bulk material. 
To test these ideas we suggest that two particle position 
correlations should be measured at smaller particle sep- 
arations where a oo l/r. Correlations should then be 
sensitive to the particle's local environment. 
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